On a Proper Definition of Spin Current 
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The conventional definition of spin current is incomplete and unphysical in describing spin trans- 
port in systems with spin-orbit coupling. A proper and measurable spin current is established in 
this study, which fits well into the standard framework of near-equilibrium transport theory and has 
the desirable property to vanish in insulators with localized orbitals. Experimental implications of 
our theory are discussed. 
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A central theme of spintronics research is on how to 
generate and manipulate spin current as well as to exploit 
its various effects [jj, |2j . In the ideal situation where spin 
(or its projection along a direction) is conserved, spin 
current is simply defined as the difference between the 
currents of electrons in the two spin states. This concept 
has served well in early study of spin-dependent transport 
effects in metals. The ubiquitous presence of spin-orbit 
coupling inevitably makes the spin non-conserved, but 
this inconvenience is usually put off by focusing one's 
attention within the so-called spin relaxation time. In 
recent years, it has been found that one can make very 
good use of spin-orbit coupling, realizing electric control 
of spin generation and transport dHH 0, IE • The 
question of how to define the spin current properly in the 
general situation therefore becomes urgent. 

In most previous studies of bulk spin transport, it has 
been conventional to define the spin current simply as 
the expectation value of the product of spin and veloc- 
ity observables. Unfortunately, no viable measurement 
is known to be possible for this spin current. The recent 
spin-accumulation experiments [a, |2| do not directly de- 
termine it, and there is no deterministic relation between 
this spin current and the boundary spin accumulation, as 
demonstrated in Fig. ^ 

In fact, the conventional definition of spin current suf- 
fers three critical flaws that prevent it from being rel- 
evant to spin transport. First, this spin current is not 
conserved. This issue alone has motivated a number of 
alternative definitions recently [1(1 Ull Il2j|. Second, this 
spin current can even be finite in insulators with local- 
ized eigenstates only, so it cannot really describe trans- 
port [13 ■ Finally, there does not exist a mechanical or 
thermodynamic force in conjugation with this current, 
so it cannot be fitted into the standard near-equilibrium 
transport theory One consequence is that one cannot es- 
tablish an Onsager relation linking the spin current with 
other transport phenomena. 

In this Letter, we try to establish a proper definition 
of spin current free from all the above difficulties, which 
is found to be possible for systems where spin generation 
in the bulk is absent due to symmetry reasons. Our new 
spin current is given by the time derivative of the spin 
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Figure 1: An example demonstrating the irrelevance of the 
conventional spin current to spin transport. Left: A macro- 
scopic system consisting of a uniform distribution of micro- 
scopic boxes (gray rectangles) in which an electron is con- 
fined. Right: The inside structure of the box, where the walls 
flips the spin of electron upon collision. When the system is 
driven to such a state schematized at the right, the macro- 
scopic average of the conventional spin current is nonzero. 
However, because all electrons are localized, no spin trans- 
port or boundary spin accumulation can occur. Tb denotes 
the total boundary spin torque due to the spin-flip scatter- 
ing. It is easy to show the corresponding spin torque dipole 
P T = —TbIb exactly cancels the conventional spin current. 
As a result, the spin current J~ s defined in Eq. J1J is zero. 



displacement (product of spin and position observables), 
which differs from the conventional definition by a torque 
dipole term. The torque dipole term is first found in a 
semiclassical theory | l 1| , whose impact on spin transport 
has been further analyzed to assess the importance of the 
intrinsic spin Hall effect |l2j |. In this work, we provide a 
quantum mechanical description of this term within the 
linear response theory for transport. Apart from showing 
its consequence in making the spin current conserved, 
we also reveal two additional properties: The new spin 
current vanishes identically in insulators with localized 
orbitals, and is in conjugation with a force given by the 
gradient of the Zeeman field or spin-dependent chemical 
potential. Together with conservation, these properties 
are crucial to establish the new spin current as the proper 
description for spin transport, and they also provide a 
firm foundation for various methods for its measurement. 

Based on general quantum mechanical principle, one 
can derive a continuity equation relating the spin, current 
and torque densities as follows, 



dS* 

at 



T, 



(1) 



2 



The spin density for a particle in a (spinor) state ip(r) 
is denned by 5 z (r) = tp' (r)s z ip(r), where s z is the spin 
operator for a particular component (z here, to be spe- 
cific). The spin current density here is given by the con- 
ventional definition J s (r) — Rei/'^(r)|{v, s z }ip(r), where 
v is the velocity operator, and {, } denotes the anticom- 
mutator. The right hand side of the continuity equation 
is the torque density defined by T z (r) = Re^ (r)r^p(r), 
where f = ds z /dt = (l/ih)[s z , H], and H is the Hamil- 
tonian of the system. These definitions can be easily 
restated in a many-body language by regarding the wave 
functions as field operators and by taking expectation 
value in the quantum state of the system. The presence 
of the torque density T z reflects the fact that spin is not 
conserved microscopically in systems with spin-orbit cou- 
pling. 

It often happens, due to symmetry reasons, that the 
average torque vanishes for the bulk of the system, i.e., 
(1/V) JdVT z (r) = 0. This is true to first order in the 
external electric field for any samples with inversion sym- 
metry. Also, one is often interested in a particular com- 
ponent of the spin, and the corresponding torque compo- 
nent can vanish in the bulk on average even for samples 
without the inversion symmetry. This is certainly true 
for the many models used for the study of spin Hall ef- 
fect |a, Li llM , and for the experimental systems used to 



detect the effect so far [8|,|9|- For such systems, where 
the average spin torque density vanishes in the bulk, we 
can write the torque density as a divergence of a torque 
dipole density, 

T z (r) = -V-P T (r). (2) 
Moving it to the left hand side of JIJ, we have 

dS. 
~dt 
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which is in the form of the standard sourceless continuity 
equation. This shows that the spin in conserved on av- 
erage in such systems, and the corresponding transport 
current is: 
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We note that there is still an arbitrariness in defin- 
ing the effective spin current because Eq. does not 
uniquely determine the torque dipole density P T from the 
corresponding torque density T z . We can eliminate this 
ambiguity by imposing the physical constraint that the 
torque dipole density is a material property that should 
vanish outside the sample. This implies in particular that 
/ dV P r = - JdVrV-P T = J dVvT z {v). It then follows 
that, upon bulk average, the effective spin current den- 
sity can be written in the form of J ' s = Retlj*(r)J s tp(r), 
where 



Js 



d(rs z ) 
dt 



is the effective spin current operator. Compared to the 
conventional spin current operator, it has an extra term 
r(ds z /dt), which accounts the contribution from the spin 
torque. 

Because the new spin current is given as a time deriva- 
tive, it must vanish in an eigen energy state in which the 
spin displacement operator is well defined, which is the 
case if the state is localized. Elementary perturbation 
theory shows that the spin current vanishes in such a 
system even in the presence of a weak electric field. In- 
deed, for spatially localized eigenstates, we can evaluate 
the spin transport coefficient as, 

lm{e\d(rs z )/dt\£'){£'\v\i) 
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where fi is the equilibrium occupation number in the ^-th 
state. Here, we have used (£|d(rs z )/dt|^') = (—i/h)(ec — 
e t ){l\rs z \£') and (£'\v\£) = (-i/H){et - ee)(t\i\t). The 
involved matrix elements are all well defined between spa- 
tially localized eigenstates. 

Defined as a time derivative of the spin displacement 
operator is z , the new spin current has a natural conju- 
gate force F s , the gradient of the Zeeman field or of a 
spin dependent chemical potential, which can be mod- 
eled as an external perturbation V = — F s • (rs z ) [lif . 
The energy dissipation rate for the spin transport can be 
written as dQ/dt = J s ■ F s . It immediately suggests a 
thermodynamic way to determine the spin current by si- 
multaneously measuring the Zeeman field gradient (spin 
force) and the heat generation. 

Moreover, Onsager relations can now be established. 
For example, in the presence of both electric and spin 
forces, the linear response of spin and charge currents 
may be written in the following manner, 
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where J s is spin current and J c denotes charge current. 
a ss and cr cc are the spin-spin and charge-charge conduc- 
tivity tensors respectively. The off diagonal block a sc 
denotes spin current response to an electric field (spin 
Hall effect), and a cs denotes charge current response to 
a spin force (inverse spin- Hall effect) [l7|. The Onsager 
reciprocity dictates a general relation between the off- 
diagonal blocks (assuming time reversal symmetry): 
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where the extra minus sign originates from the odd time- 
reversal parity of the spin displacement operator rs z |l8| . 
We note that the Onsager relation Eq. JHJ can only be 
established when the currents are defined in terms of the 
time derivative of displacement operators conjugate to 



3 



the forces. In the case of spin force driving, we have 
the spin displacement operator s z r, whose time deriva- 
tive corresponds to our new spin current Eq. (J5J , not the 
conventional spin current. This Onsager relation can also 
be directly verified using the linear response theory. 

The existence of the Onsager relations makes the 
electric measurement of the spin current viable. Previous 
theoretical proposal had suggested that the spin current 
can be determined by measuring the transverse voltage 
generated by a s pin current passing through a spin 
Hall device 0, lil llflj. However, to do that, the spin 
Hall coefficient cr*° of the measuring device must be 
known in prior because the spin current is determined 
from J sx — (j s x c y E y . With the Onsager relation, a^y 
can be derived from the corresponding inverse spin 
Hall coefficient, and the latter can be determined by a 
measurement of the charge current and the Zeeman field 
gradient. 

After these general considerations, we now show how 
to evaluate the spin Hall conductivity based on the new 
definition of spin current. The torque dipole density can 
be determined unambiguously as a bulk property within 
the theoretical framework of linear response. Consider 
the torque response to an electric field at finite wave vec- 
tor q, T z (q) = x(q) • E(q). Based on Eq. © which 
implies T z {c\) = — iq • P r (q), we can uniquely determine 
the dc response (i.e., q — > 0) of the spin torque dipole: 

P r = Rc{zV q [x(q)-E]} q=0 . (9) 

Here we have utilized the condition x(0) = 0, i.e., there is 
no bulk spin generation by the electric field. Combining 
Eq. (@J and ©, we can then determine the electric-spin 
transport coefficients for the new definition of spin cur- 
rent: 

o%, = < + (10) 

where a s ® v is the conventional spin-transport coefficient 
that is the focus of most of previous studies, and 

^ v = R.e[j^x^(q)] q =o , (if) 

is the contribution from the spin torque dipole. Standard 
Green function or many-body techniques can be used to 
evaluate this linear response for systems with arbitrary 
disorder and interactions between the carriers. 

The spin-Hall coefficients for a few semiconductor 
models including effects of disorder are now considered by 
Sugimoto et al. based on our new definition of spin cur- 
rent [2(|, who found results dramatically different from 
the conventional spin Hall conductivities. For example, it 
is found that the spin Hall conductivity depends explic- 
itly on the scattering potentials for the two dimensional 
Rashba models with k-linear or k-cubic spin-orbit cou- 
pling. For the k-cubic model, the conventional spin Hall 



conductivity is robust against disorder, but this is not 
so if the new spin current definition is adopted. It then 
implies that at smooth boundaries, spin accumulation in 
such systems are of extrinsic nature. 

The conservation of our new spin current allows one to 
consider spin transport in the bulk without the need of 
laboring explicitly a spin torque (dipole density) which 
may be generated by the electric field. One can think 
of spin transport for systems with strong spin-orbit cou- 
pling in the "usual" sense established for systems with 
weak spin-orbit coupling. For example, it has been cus- 
tomary to link spin density and spin current through the 
following phenomenological equation of spin continuity, 

+ = (12) 

at t s 

where r s is the spin relaxation time, and the spin current 
has the form J s = aE — D S VS Z . This makes sense 
only if our new spin current is used in the calculation 
of spin Hall conductivity a, otherwise an extra term of 
field-generated spin torque must be added [23j. 

Equation l|12|) can serve as the basis to determine 
the spin accumulation at a sample boundary, which is 
of much current interest. Consider a system having a 
smooth boundary produced by a slowly varying confin- 
ing potential. We assume that the length scale of varia- 
tion is much larger than the mean free path, so that the 
above continuity equation may be applied locally. By in- 
tegrating from the interior to the outside of the sample 
boundary, we obtain a spin accumulation per area with 
S z = J^ ulk T s 0], where r s is the spin relaxation time. 
We emphasize that the transport spin current responsi- 
ble for the boundary spin accumulation should be JT S 
instead of the conventional spin current J s . 

For sharp boundaries, the spin continuity equation 
alone cannot yield a unique relationship between the 
spin current from the bulk and spin accumulation at the 
boundary. For a perfectly reflecting sharp wall in the 
case of strong spin-orbit coupling, the boundary spin ac- 
cumulation seems to be determined by the conventional 
spin current from the bulk |2l|. However, such a relation- 
ship is altered for other boundary conditions [2^]. This 
calls for well controlled experiments to explicitly elimi- 
nate the influence from the boundary condition. On the 
other hand, for the generic class of smooth boundaries 
discussed above, there is a unique relationship between 
spin accumulation and the spin current, provided one 
uses our new definition. 

The real advantage of our new definition of spin current 
lies in the fact that it provides a satisfactory description 
of spin transport in the bulk. With our new spin current, 
one can now use the spin continuity equation (12) to dis- 
cuss spin accumulation in the bulk, e.g., by generating 
a non-uniform electric field or spatially modulating the 
spin Hall conductivity. Our new spin current vanishes in 
Anderson insulators either in equilibrium or in a weak 
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electric field, which enables us to predict zero spin accu- 
mulation in such systems. More importantly, it posses a 
conjugate force (spin force), so that spin transport can 
be fitted into the standard formalism of near equilibrium 
transport. The conventional spin current does not have a 
conjugate force, so it makes no sense even to talk about 
energy dissipation from that current. The existence of 
a conjugate force is crucial for the establishment of On- 
sager relations between spin transport and other trans- 
port phenomena, and its measurement will be important 
to thermodynamic and electric determination of the spin 
current. 

We gratefully acknowledge discussions with S. Zhang, 
Z. Yang, D. Culcer, J. Sinova, S.-Q. Shen and A.H. Mac- 
Donald. QN and DX were supported by DOE (DE-FG03- 
02ER45958) and the Welch Foundation, and JRS and 
PZ were supported by NSF (DMR-0306239 and DMR 
0404252). 



* These authors contributed equally to this work. 
[1] G.A. Prinz, Science 282, 1660 (1998); S.A. Wolf et al, 

Science 294, 1488 (2001). 
[2] I. Zutic, J. Fabian, and S.D. Sarma, Rev. Mod. Phys. 76, 

323 (2004). 

[3] M.I. Dyakonov and V.I. Perel, JETP 33, 1053 (1971). 
[4] J. E. Hirsch, Phys. Rev. Lett. 83, 1834 (1999). 
[5] S. Zhang, Phys. Rev. Lett. 85, 393 (2000). 
[6] S. Murakami, N. Nagaosa, and S.C. Zhang, Science 301, 
1348 (2003). 

[7] J. Sinova et al., Phys. Rev. Lett. 92, 126603 (2004). 
[8] Y.K. Kato et al, Science 306, 1910 (2004). 
[9] J. Wunderlich, B. Kaestner, J. Sinova, T. Jungwirth, 
Phys. Rev. Lett. 94, 047204 (2005). 



[10] P.-Q. Jin, Y.-Q. Li and F.-C. Zhang, preprint, 
cond-mat/0502231 

[11] S. Murakami, N. Nagaosa, and S.-C. Zhang, Phys. Rev. 
B 69, 235206 (2004). 

[12] S. Zhang and Z. Yang, Phys. Rev. Lett. 94, 066602 
(2005). In reaching their conclusion Eq. (17), they made 
an implicit assumption that (k\\[rs, H]\k\) = 0. How- 
ever, because fcA) is an extended state, this relation is in 
general not valid. We further note that {k\\[r, H]\k\) ^ 
is why we can see non-vanishing conductivity in metals. 

[13] E.I. Rashba, Phys. Rev. B 68, 241315(R) (2003). Ac- 
tually (v x o y - Vy<Jx) = -(h/m)[dE(ot R ) / da R ] - 2a R is 
non-zero even in the presence of a confinement potential. 

[14] D. Culcer et al, Phys. Rev. Lett. 93, 046602 (2004). 

[15] J. Schliemann and D. Loss, Phys. Rev. B 69, 165315 
(2004); 71, 085308 (2005). 

[16] I. Zutic, J. Fabian, S. Das Sarma, Phys. Rev. Lett. 88, 
066603 (2002); J. Fabian, I. Zutic, S. Das Sarma, Phys. 
Rev. B 66, 165301 (2002). 

[17] P. Zhang and Q. Niu, unpublished, cond-mat/0406436 

[18] H.B. Casimir, Rev. Mod. Phys. 17, 843 (1945). 

[19] E.M. Hankiewicz, L.W. Molenkamp, T. Jungwirth, J. 
Sinova, Phys. Rev. B 70, 24130RR) (2004). 

[20] N. Sugimoto, S. Onoda, S. Murakami and N. Nagaosa, 
preprint, cond-mat/0503475 The calculation of the spin 
Hall coefficients for a number of semiconductor mod- 
els in the clean limit was first carried out by us (see 
cond-mat/0503505 v. 4). The dramatically different re- 
sults were also found. 

[21] K. Nomura, J. Wunderlich, Jairo Sinova, B. Kaest- 
ner, A.H. MacDonald, T. Jungwirth, preprint, 
|cond-mat/ 0508532 

[22] Wang-Kong Tse, J. Fabian, I. Zutic, S. Das Sarma, 
preprint, cond-mat /0508076 

[23] In Ref . fl2| , the extra contribution of field generated spin 
torque was also considered. However, they did not reach 
the correct spin current definition due to an invalid as- 
sumption mentioned in the remark of p^|. 



